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ABSTRACT 


A  unidirectional  fiberglass  epoxy  composite  specimen  is 
modelled  as  a  homogeneous  transversely  isotropic  continuum  plate 
medium.  Acousto-ultrasonic  non-contact  input-output  characteri¬ 
zation  by  tracing  SV  waves  in  the  continuum  is  studied  theoreti¬ 
cally  with  a  transmitting  and  a  receiving  transducer  located  on 
the  same  face  of  the  plate.  The  isotropic  plane  of  the  equi¬ 
valent  continuum  plate  model  lies  in  the  midplane  of  the  plate 
and  is  parallel  to  the  top  and  bottom  faces  of  the  plate. 


The  single  reflection  problem  for  an  incident  SV  wave  at  a 
stress -free  plane  boundary  in  a  semi -infinite  transversely  iso¬ 
tropic  medium  whose  isotropic  plane  is  parallel  to  the  plane 
boundary  is  analyzed  first.  For  all  angles  of  incidence,  the 
angle  of  reflection  of  the  SV  wave  is  equal  to  the  angle  of  inci¬ 
dence  of  the  SV  wave.  It  is  found  that  an  obliquely  incident  SV 
wave  results  in  a  reflected  SV  wave  and  a  reflected  P  wave  for  an 


angle  of  incidence  of  the  incident  SV  wave  less  than  the  critical 
angle  52®.  when  the  angle  of  incidence  of  an  incident  SV  wave  is 
equal  to  or  greater  than  the  critical  angle,  there  exists  only  an 


SV  wave  in  the  medium  as  the  reflected  P  wave  degenerates  into  a 
surface  wave  travelling  parallel  to  the  plane  boundary.  The 
amplitude  of  the  surface  wave  decays  exponentially  with  the  per¬ 
pendicular  distance  from  the  plane  boundary.  The  amplitude 


*Qn  For 

.•RA&I 

8 

□  m 

.?ced 

□  1 

By 

Avail  ■■ 

Blst  Specif 

A  ^  1  i  I  ' 


ratios  of  the  reflected  P  and  SV  waves  to  the  incident  SV  wave  as 
a  function  of  the  angle  of  incidence  are  plotted  from  zero  to  the 
critical  angle.  The  amplitude  ratio  of  the  reflected  SV  wave  is 
found  to  be  minus  one  when  the  angle  of  incidence  is  equal  to  or 
greater  than  the  critical  angle.  The  balance  in  energy  flux  nor¬ 
mal  to  the  plane  boundary  is  checked. 

Accordingly,  the  delay  time  for  wave  propagation  between  the 
transmitting  and  the  receiving  transducers  is  computed  as  if  the 
SV  waves  were  propagating  in  a  half  space.  It  is  found  that  the 
directional  dependence  of  the  phase  velocity  of  the  SV  wave  pro¬ 
pagating  in  the  transversely  isotropic  medium  has  a  significant 
effect  on  the  delay  time,  as  opposed  to  the  directional  indepen¬ 
dence  of  the  phase  velocity  of  a  shear  wave  propagating  in  an 
isotropic  medium. 

The  displacements  associated  with  the  SV  wave  in  the  plate 
and  which  may  be  detected  by  the  non- contact  receiving  transduc¬ 
er  are  approximated  by  an  asymptotic  solution  for  an  infinite 
transversely  isotropic  medium  subjected  to  a  harmonic  point  load. 
The  polar  diagrams  for  the  directivity  function  of  the  stresses 
due  to  SV  waves  in  the  plate  are  shown  at  frequencies  of  0.75, 
1.50  and  2.25  MHz. 
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This  study  enhances  the  quantitative  understanding  of  acous- 
to-ultrasonic  non-destructive  evaluation  (NDE)  parameters  such  as 
the  stress  wave  factor  (SWF)  and  wave  propagation  in  fiber  rein¬ 
forced  polymeric,  ceramic  or  metallic  composites,  which  can  be 
modelled  as  transversely  isotropic  media. 
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INTRODUCTION 


Fiber  reinforced  composite  materials  are  attractive  materials 
for  aerospace  applications  because  of  their  high  specific 
mechanical  properties.  It  has  been  shown  that  many  composites, 
such  as  fiberglass  epoxy  composites  or  fiber  reinforced  ceramics 
as  shown  in  Fig,  1,  may  be  modelled  as  a  homogeneous  transversely 
isotropic  continuum  [1].  In  this  work,  acousto-ultrasonic  (AU) 
non-contact  input-output  characterization  of  a  homogeneous  trans¬ 
versely  isotropic  elastic  plate  is  investigated  by  tracing  SV 
waves . 

First,  the  single  reflection  problem  of  an  incident  SV  wave 
at  a  stress -free  plane  boundary  in  a  semi- inf inite  transversely 
isotropic  medium  whose  isotropic  plane  is  parallel  to  the  plane 
boundary  is  considered.  At  such  boundaries,  the  conditions  for 
the  existence  of  wave  mode  conversion,  critical  angle  phenomena, 
reflected  surface  wave,  the  angles  of  reflection  of  the  reflected 
waves,  and  the  amplitude  ratios  of  the  reflected  waves  to  the 
incident  wave  are  derived. 

Second,  the  SV  wave  input -output  relations  are  derived  when 
multiple  reflections  occur  at  the  top  and  the  bottom  faces  of  the 
plate.  The  delay  time  between  input  and  output  versus  the  dis- 
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tance  separating  the  transmitting  and  receiving  transducers  is 
analyzed.  The  directivity  functions  of  the  stresses  associated 
with  the  SV  waves  are  computed.  And,  the  output  displacement  at 
the  non-contact  receiving  transducer  is  approximated  by  an 
asymptotic  solution. 

This  investigation  should  enhance  the  quantitative  understan¬ 
ding  of  AU  NDE  parameters  such  as  the  stress  wave  factor.  It 
also  provides  the  potential  for  assisting  in  the  development  of 
better  NDE  schemes  utilizing  the  SWF  as  well  as  other  AU  parame¬ 
ters  for  the  characterization  of  fiber  reinforced  polymeric, 
metal  and  ceramic  composites . 
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SINGLE  REFLECTION  PROBLEM  AT  STRESS -FREE  PLANE  BOUNDARY  IN 
SEMI -INFINITE  TRANSVERSELY  ISOTROPIC  MEDIUM  WHOSE  ISOTROPIC  PLANE 
IS  PARALLEL  TO  PLANE  BOUNDARY  FOR  INCIDENT  SV  WAVE 
1 .  Reflected  P  and  SV  Waves 


For  a  homogeneous,  linearly  elastic  transversely  isotropic 
continuum,  the  number  of  independent  elastic  constants  is  five 
[1].  Define  a  coordinate  system  (x,  y,  z)  for  a  semi- inf inite 
transversely  isotropic  medium  whose  isotropic  plane  is  parallel 
to  the  plane  boundary  where  the  reflection  occurs  as  follows:  the 
plane  boundary  contains  the  x  and  y  axes,  and  the  z  axis  is  the 
zonal  axis  of  the  medium,  which  is  in  the  direction  parallel  to 
the  fiber  direction  shown  in  Fig.  1.  Also,  see  Fig.  2.  The  gen¬ 
eralized  Hooke's  law  is  written,  relative  to  the  (x,  y,  z)  coor¬ 
dinate  system,  as  [1] 

^xx  “  Ciiu,x  +  y  +  Ci3W,z 

^yy  ""  +  C]_2^v,y  +  Cx3W,z 

^zz  ^  Ci3U,x  +  Ci3V,y  +  C33W,z 

^xz  =  C44(u,z  +  w,x) 

^yz  ^  C44(v,z  +  w,y) 

'^xy  “ 

where  r^g  (r,  s=x,  y  and  z)  are  the  normal  (r  ^  s)  and  shear 


(1) 
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(r  5^  s)  stresses  with  respect  to  the  coordinate  system  (x,  y,  z)  ; 
u,  V  and  w  are  the  displacement  components  of  a  point  in  the 
medium  along  the  x,  y  and  z  axes,  respectively;  ”  denotes  par¬ 
tial  differentiation  with  respect  to  the  variable  which  follows; 
and  ^12 »  ^13 »  ^33  ^44  five  independent  elastic 

constants  where  Cgg  =  -  ^12)- 

Let  a  progressive  wave  be  represented  as  [2] 

(u,  V,  w)  =  A  (P^,  Py,  Pz)  exp{ia)(SxX  +  Syy  +  S^z  -  t) }  (2) 

where  Sy  and  S2  are  the  components  of  the  slowness  vector, 
which  is  in  the  same  direction  as  the  normal  to  the  wavefront  and 
whose  magnitude  is  equal  to  the  reciprocal  of  the  magnitude  of 
the  phase  velocity  [1],  along  the  x,  y  and  z  axes,  respectively; 
P^,  Py  and  P2  are  the  components  of  a  unit  vector  of  particle 
displacement  along  the  x,  y  and  z  axes,  respectively;  A  is  the 
amplitude  of  particle  displacement;  t  denotes  time  and  w  denotes 
radian  frequency.  It  follows  from  Eqs.  (1)  and  (2)  that  the 
stresses  can  be  represented  as 

^xx  “  ^12Sy^y  Ci3S2P2]exp{  iwCS^x  +  Syy  +  S^z  -  t) } 

Tyy  =  iwA[Cj^2SxPx  +  +  Ci3S2Pz]exp{  iwCS^x  +  Syy  +  S^z  -  t)  } 

r^z  =  iwA[Ci3SxPx  +  Ci3SyPy  +  C33S2Pz]exp{ iwCS^x  +  Syy  +  S^z  -  t) ) 
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(3) 


Txz  -  iwAtC^^S^Pjj  +  C44SxPz]exp{iw(SxX  +  Syjr  +  S^z  -  t) ) 

^yz  “  iwA[C44S2Py  +  C44SyP2]exp{ iwCS^x  +  Syy  +  S^z  -  t) ) 

Txy  -  iwA[C65SyPx  +  C66SjjPy]exp{iw(SxX  +  Syy  +  S2Z  -  t) ) 

The  stress  boundary  conditions  on  the  stress -free  plane 
boundary  require  that  [2] 


+ 

r  (R) 

-  0 

+ 

-yz<‘'> 

■=  0 

(4) 

+ 

-  0 

where  “  x,  y  and  z)  represent  stresses  on  the  plane 

boundary  associated  with  the  incident  SV  wave,  and  (r  =>  x, 

y  and  z)  represent  stresses  on  the  plane  boundary  associated  with 
the  reflected  waves. 

In  order  to  satisfy  Eq.  (4),  it  is  required  [2]  that  the 
frequency,  w,  of  the  reflected  wave  be  equal  to  that  of  the 
incident  wave  and  that  [2] 

S,<I)  -  S,W 

SyU)  -  SyfR) 
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As  a  result  of  Eq.  (5),  the  slowness  vectors  of  the  incident  and 
reflected  waves  lie  in  a  plane  called  the  plane  of  incidence. 

This  analysis  can  be  simplified  by  assuming  that  the  plane  of 
incidence  is  the  x  *  0  plane;  that  is,  the  slowness  vectors  of 
the  incident  and  reflected  waves  are  in  the  x  0  plane,  as  shown 
in  Fig.  2.  Then  it  follows  from  Eqs.  (5)  that 

-  0  (6) 

It  has  been  shown  [3]  that,  except  along  the  principal  mate¬ 
rial  axes,  P  waves  or  SV  waves  travelling  in  a  plane  containing 
the  zonal  axis,  z  axis,  of  a  transversely  isotropic  medium  are 
quasi -longitudinal  and  quasi- transverse,  respectively;  that  is, 
in  general,  the  components  of  the  unit  vector  of  particle  dis¬ 
placement  for  either  the  P  or  SV  waves  along  the  y  and  z  axes, 

Py  and  P^,  do  not  vanish;  whereas  the  components  along  the  x 
axis,  Pj^,  do  vanish.  Therefore,  it  follows  from  Eqs.  (3)  and  (6) 
that  stresses  associated  with  P  and  SV  waves  are 

^xz  -  ’’xy  “  0 

(7) 

^xx  ^  ^yy  ^  ^zz  ^  ^yz  ® 

It  has  also  been  shown  that  an  SH  wave  travelling  in  a  plane 
containing  the  zonal  axis.z  axis,  of  a  transversely  isotropic 
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medium  possesses  a  transverse  displacement  only,  that  is,  for  the 
coordinates  in  Fig.  2,  (P^,  Py,  Pz)  -  (1,  0,  0)  [3].  Therefore, 
it  follows  from  Eqs.  (3)  and  (6)  that  the  stresses  associated 
with  SH  waves  are 


^yz  ^xx  ^3ry  ~  ^zz  ~  ® 


Txy  ^  0;  Tjjz  ^  0 


Assume  that  an  SV  wave  is  incident  on  the  plane  boundary,  the 
x-y  plane  in  Fig.  2.  It  follows  from  Eqs.  (4)  and  (7)  that 


^xz('>  =  0 


As  a  result  of  Eq.  (9),  it  is  known  from  Eq.  (4)  that  is 

equal  to  zero.  This  means  that  no  SH  wave  will  be  reflected  back 
into  the  medium  because  a  reflected  wave  of  the  SH  type  would 
result  in  a  nonzero  value  of  the  stress  r^^z  •  So,  Eqs.  (4) 
reduce  to 


^zz^'>  -  rzz(^>  =  0 
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Since  either  a  reflected  P  wave  or  a  reflected  SV  wave  results  in 
nonzero  values  of  the  stresses  and  therefore 

concluded  from  Eqs .  (7)  and  (10)  that  both  a  P  wave  and  an  SV 
wave  may  be  reflected  back  into  the  medium. 

2 .  Slowness  Surface  for  P  and  SV  Waves 

The  equations  of  motion  relative  to  the  coordinate  system  (x, 
y,  z)  are  [1] 

^xx,x  ^xy,y  ^xz,z  ""  pu,tt 

’’xy.x  +  ’■yy.y  ’"yz.z  =  P^.tt  (11) 

’’xz.x  ^yz,y  ’’zz.z  ”  pw,tt 

where  the  body  forces  are  identically  zero  for  the  homogeneous 
solution. 

It  follows  from  Eqs.  (1),  (2)  and  (11)  that  the  following 
equations  of  motion  are  obtained: 

t^ll^x^  ^66^y^  ^44^z^  “  ^^12  ^66^^x^y^y 

+  (Ci2  +  “  0 
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(12) 


+  (C]^3  4-  C^^)SyS^?^  -  0 

(C]^3  +  ^  ^^13  ■*■  ^44)SySz^y  [C44(Sx^  ^  ^y^) 

^33^z^  "PI^z  *  ^ 

The  condition  for  the  existence  of  the  plane  wave  solution  is 
expressed  by  setting  the  determinant  of  the  matrix  of  the  coef¬ 
ficients  of  Py  and  P2  in  Eq.  (12)  equal  to  zero  [1]: 

CCiiS^2  ^  c^Sy2  *  C44S,2  -p]  "  C66)SxSy  (C13  +  044)8^$, 

<Ci2  +  +  c,,Sy2  +  0448^2  .^]  (Ci3  +  C44)8y82 

^^^13  *  ^*^13  *  *  ^3Zh^ 

(13) 

By  expanding  Eq.  (13),  three  sheets  of  slowness  surface  are 
obtained.  The  slowness  surface  for  an  SV  wave  is  given  in  [3]  as 

Cii  +  C44  C44  +  C33 

( - )(Sx2  +  S  2)  +  ( - )  5^2 

2  ^  2 

r  r 

-  1/2-  (Cii  -  C44)(Sx2  +  S  2)  +  (C33  - 

[[  J  (1^) 

+  4(Sx2  +  SyhSz'^  (Cii  -  C44)(C33  -  C44) 


-  (Ci3  +  C44)2 


-  P 


where  S^,  Sy  and  are  the  components  of  the  slowness  vector  for 
an  SV  wave  along  the  x,  y  and  z  axes,  respectively;  p  is  density; 
and  Cii,  ^13 y  C33  and  C44  are  elastic  constants.  Similarly,  the 
slowness  surface  for  a  P  wave  is  [3] 


+  l/2||^(Cii  -  C44)(S^2  +  Sy2)  +(C33  -  €44)8^2 

+4(8^2  +  83,2)8^2  -  C44)(C33  +  C44) 

-  (Ci3  +  C44)2j|  ^  =  p 


where  Sy  and  S^  are  the  components  of  the  slowness  vector  for 
a  P  wave  along  the  x,  y  and  z  axes,  respectively. 


One  quadrant  of  the  intersection  of  the  slowness  surface  of 
an  SV  wave  in  the  unidirectional  fiberglass  epoxy  composite  shown 
in  Fig.  1  with  the  plane  x  *  0  and  one  quadrant  of  the  intersec¬ 
tion  of  the  slowness  surface  of  a  P  wave  in  the  unidirectional 
fiberglass  epoxy  composite  with  the  plane  x  =«  0  are  shown  from 
Eqs.  (14)  and  (15)  in  Fig.  3;  where  the  numerical  values  of  the 
elastic  constants  and  density  given  in  [1]  for  the  unidirectional 
fiberglass  epoxy  composite  are  used  and  are  as  follow:  = 

10.581  X  10%/m2,  C;l3  -  4.679  x  10%/m2,  C33  -  40.741  x  109N/m2, 
C44  =  4.422  X  10^N/ra2,  and  p  -  1850  kg/m^ . 
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^ •  Angle  of  Reflection 


It  follows  from  Eqs.  (5)  and  (6)  that  the  y-component  of  the 
slowness  vector  of  an  incident  SV  wave  is  equal  to  the  y- compo¬ 
nent  of  the  reflected  SV  wave  as  well  as  the  y-component  of  the 
reflected  P  wave.  Accordingly,  the  relation  between  the  y-com¬ 
ponent  of  the  slowness  vector  of  an  incident  SV  wave  and  the 
y-component  of  the  slowness  vectors  of  the  reflected  SV  and  P 
waves  is  given  as 

Syd)  -  Sy(SV)  .  Sy<P)  -  b  (16) 

where  represents  the  y-component  of  the  slowness  vector  of 

an  incident  SV  wave;  Sy^^^)  represents  the  y-component  of  the 
slowness  vector  of  the  reflected  SV  wave;  Sy^^^  represents  the 
y-component  of  the  slowness  vector  of  the  reflected  P  wave;  and  b 
is  a  common  constant,  as  shown  in  Fig.  3. 

It  follows  from  Eqs.  (6),  (14)  and  (16)  that  the  relation 
between  the  z- component  of  the  slowness  vector  of  an  incident  SV 
wave  and  that  of  the  reflected  SV  wave  is 

s^(l)  .  .  s^(SV)  (17, 

The  minus  sign  is  due  to  the  fact  that  the  slowness  vector  of  an 
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incident  SV  wave  points  out  of  the  medium,  whereas  the  slowness 
vector  of  the  reflected  SV  wave  points  into  the  medium,  as  shown 
in  Fig.  2.  Consequently,  the  value  of  the  z-component  of  the 
slowness  vector  of  an  incident  SV  wave,  S^^^^  is  negative, 
whereas  that  of  the  reflected  SV  wave,  is  positive. 

The  angle  of  reflection  is  defined  as  the  angle  between  the 
slowness  vector  of  a  reflected  wave,  either  type  P  or  SV,  and  the 
normal  to  the  plane  boundary  where  the  reflection  occurs.  Simi¬ 
larly,  the  angle  of  incidence  is  defined  as  the  angle  between  the 
slowness  vector  of  an  incident  SV  wave  and  the  normal  to  the 
plane  boundary,  as  shown  in  Fig.  2.  Therefore,  the  angle  of 
reflection  of  a  reflected  SV  wave  is  defined  as 

«  tan“^(Sy  (SV)/s^(SV))  (18) 

and  the  angle  of  incidence  of  an  incident  SV  wave  is  defined 
as 

9i  =  tan-l(Sy(I)/-S2(I))  (19) 

It  follows  from  Eqs.  (16)  through  (19)  that,  for  all  angles,  the 
angle  of  incidence  of  an  incident  SV  wave  is  equal  to  the  angle 
of  reflection  of  the  reflected  SV  wave,  as  shown  in  Fig.  2. 
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However,  the  angle  of  reflection  of  the  reflected  P  wave  is 
not  equal  to  the  angle  of  incidence  of  the  incident  SV  wave.  For 
^  given  value  of  b  in  Eq.  (16),  two  values  of  the  z-component  of 
the  slowness  vector,  S^,  of  equal  magnitude  but  opposite  sign 
(for  a  P  wave  travelling  in  the  plane  x  -  0  in  the  transversely 
isotropic  medium)  can  be  obtained  from  the  slowness  surface  for 
the  P  wave  by  substituting  Eqs.  (6)  and  (16)  into  Eq.  (15).  The 
positive  z-component  of  the  slowness  vector  corresponds  to  the 
reflected  P  wave  and  is  denoted  as  Similarly,  for  a  given 

vQ-lue  of  b,  there  exists  a  positive  z-component  of  the  slowness 
vector,  corresponding  to  the  reflected  SV  wave.  Accord¬ 

ingly,  for  any  given  value  of  b  in  Eq.  (16),  there  exist  a  posi¬ 
tive  z-component  of  the  slowness  vector  for  the  reflected  SV 
wave,  ^  ^ ,  and  a  positive  z-component  of  the  slowness  vector 

for  the  reflected  P  wave,  62^^).  In  fact,  the  z-components  of 
the  slowness  vectors  of  the  reflected  SV  and  P  waves,  S2^^^^ 

incident  SV  wave  which  determines  the  value  of  b  in 
Eq .  (16)  and  which  travels  in  the  plane  x  “■  0  in  the  transversely 
isotropic  continuum,  can  be  obtained  from  the  lengths  of  the  per¬ 
pendicular  lines  between  the  abscissa  representing  the  value  of  b 
and  the  intersections  with  the  two  sheets  of  the  slowness  sur¬ 
faces  for  the  reflected  SV  and  the  reflected  P  waves,  respec¬ 
tively,  as  shown  in  Fig.  3.  It  is  apparent  from  Fig.  3  that  the 
z-component  of  the  slowness  vector  of  the  reflected  SV  wave. 
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is  greater  than  that  of  the  reflected  P  wave,  for 

an  incident  SV  wave  travelling  in  the  plane  x  =  0  in  the  unidi¬ 
rectional  fiberglass  epoxy  composite  [1].  Consequently,  the 
angle  of  reflection  of  the  reflected  P  wave  9-p,  defined  similarly 
to  Eq.  (18)  as 

=  tan-l(Sy(P)/S2<P)),  (20) 

is  larger  than  the  angle  of  reflection  of  the  reflected  SV  wave 
for  an  incident  SV  wave  in  the  unidirectional  fiberglass  epoxy 
composite,  from  Eqs .  (16),  (18)  and  (20). 

The  critical  angle  is  defined  as  the  angle  of  incidence  of  an 
incident  SV  wave  at  which  the  slowness  vector  of  the  reflected  P 
wave  becomes  tangent  to  the  plane  boundary  where  the  reflection 
occurs.  For  the  specific  case  of  an  incident  SV  wave  in  the  uni¬ 
directional  fiberglass  epoxy  composite  in  Fig.  1,  it  has  been 
shown  above  that  the  angle  of  reflection  of  the  reflected  P  wave 
is  greater  than  that  of  the  reflected  SV  wave.  Thus,  the  criti¬ 
cal  angle  phenomenon  occurs  for  this  composite.  The  critical 
angle  for  an  incident  SV  wave  in  the  unidirectional  fiberglass 
epoxy  composite  shown  in  Fig.  1  is  equal  to  52°.  The  angle  of 
reflection  of  the  reflected  P  wave  versus  the  angle  of  incidence 
of  an  incident  SV  wave  in  the  unidirectional  fiberglass  epoxy 
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composite  [1]  is  plotted  from  zero  to  the  critical  angle  52®,  as 
shown  in  Fig.  4. 

4.  Amplitude  Ratios  of  Reflected  Waves  to  Incident  Wave  for 

Angle  of  Incidence  Less  Than  Critical  Angle 

It  has  been  stated  that  an  SV  wave  and  a  P  wave  will  be 
reflected  back  into  the  semi -infinite  transversely  isotropic 
medium  if  the  angle  of  incidence  of  the  SV  wave  incident  on  the 
plane  boundary  is  less  than  the  critical  angle.  The  stress 
boundary  conditions  for  an  incident  SV  wave  whose  slowness  vec¬ 
tor  is  in  the  plane  x  =*  0  can  be  rewritten,  from  Eqs.  (4)  and 
(7),  as 


,„(I)  +  .  0 


(21) 


where  and  represent  the  shear  and  normal  stresses 

associated  with  the  incident  SV  wave;  and 

represent  the  shear  and  normal  stresses  associated  with  the 
reflected  SV  wave;  and  Ty^^^^  and  represent  the  shear  and 

the  normal  stresses  associated  with  the  reflected  P  wave. 

The  shear  stress,  associated  with  an  SV  wave  of  unit 

amplitude  travelling  in  the  plane  x  =  0  and  incident  on  the  plane 
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boundary  at  the  origin,  as  shown  in  Fig.  2,  can  be  obtained  from 
Eqs .  (3)  and  (6)  as 

-  ia,(C55S2<I>Py<^>  +  CssSyd)?^^^))  (22) 

where  and  are  the  components  of  the  unit  vector  of 

particle  displacement  of  the  incident  SV  wave  along  the  y  and  z 
axes,  respectively.  Similarly,  the  normal  stress, 
associated  with  the  incident  SV  wave  of  unit  amplitude  can  be 
expressed,  from  Eqs.  (3)  and  (6)  as 

=  ic.(Ci3Sy(I)Py(I)  +  CajS^CDp^d))  (23) 

The  shear  stress,  and  the  normal  stress, 

associated  with  the  reflected  SV  wave  on  the  plane  boundary  at 
the  origin  are,  from  Eqs.  (3)  and  (6) 

ry^(SV)  ^  i<,A(SV)(C55S^(SV)py(SV)  +  C55Sy(SV)p^(SV) )  (24) 

and 

r„(SV)  .  l^(SV)(Ci3Sy(SV)pj,(SV)  ^  C33S^<SV)p^(SV) ,  (25) 

where  Py^^^^  and  P2^^^^  are  the  components  of  the  unit  vector  of 
particle  displacement  of  the  reflected  SV  wave  along  the  y  and  z 


axes,  respectively;  and  Is  the  amplitude  of  the  reflected 

SV  wave . 

The  shear  stress,  and  the  normal  stress, 

associated  with  the  reflected  P  wave  on  the  plane  boundary  at  the 
origin  are,  from  Eqs.  (3)  and  (6), 

=  iu>A(P)(C55S2(^)Py(^)  +  C55Sy(P)p2<P>)  (26) 

and 

=  iwA(P)(Ci3Sy(P)Py(P)  +  C33S2(^>Pz(Pb  (27) 

where  Py^^^  and  Pz^^^  are  the  components  of  the  unit  vector  of 
particle  displacement  of  the  reflected  P  wave  along  the  y  and  z 
axes,  respectively;  and  A^^^  is  the  amplitude  of  the  reflected  P 
wave . 

Upon  substitution  of  Eqs.  (22)  through  (27)  into  Eq.  (21), 
the  stress  boundary  conditions  for  an  incident  SV  wave  of  unit 
amplitude  travelling  in  the  plane  x  -  0  can  be  expressed  as 
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(28) 


S^(I)Py(I)  +  SyfDp^d)  +A<J’)(S2(f)Py<'')  +  Sy(P)P2<P)) 

+  A<S'')(Sj(SV)Py(SV)  +  Sy(S'')P2(SV))  -  0 

(Ci3Sy(I)Py(I)  +  C33S^(I)p^(I))  +  A<P)(Ci3Sy(r)Py(P) 

+  C33S^<P)p^(P))  +  A(SV)(C33Sy(SV)Py<SV) 

+  C33Sj(S'')p^(SV),  .  0 

The  components  of  a  unit  vector  of  particle  displacement 
along  the  y  and  z  axes  of  an  SV  wave  travelling  in  any  plane 
containing  the  zonal  axis  of  a  transversely  isotropic  medium  are 
given  in  [3],  and  when  applied  to  the  present  case  of  the 
reflected  SV  wave  travelling  in  the  plane  x  =  0,  can  be  expressed 
as 


1  /([  "1 2 

■  (^33  “  ^SV  ■  (^33  '  ^44^^z^^^^ 

^  [(Ci3  *  C44)Sy<S'')s,(SV) 

(29) 

1  /  r  r  1 2 

p^(SV)  .  (Cj3  *  C44)Sy<^'')Sj,<P''>  /•  Hjv  -  (C33  ■  C44)Sj(®''> 

^  [(=13  *  C44)Sy(='')s,(SV) 

J  J 
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where  Hctt  is  defined  as 


HSV  =  ](Cll  -  C44)Sy(SV)^  +  (C33  -  C44)Sy(SV)^ 


(Cii  -  C44)Sy<SV)^  (C33  -  C44)S2<SV)^  ^ 

-  4S  (SV)2s^(SV)2  .  C44)(C33  -  C44) 


-  (Ci3  +  044)' 


2II/2 


Similarly,  the  components  of  a  unit  vector  of  particle  displace¬ 
ment  along  the  y  and  z  axes  for  the  reflected  P  wave  travelling 
in  the  plane  x  =  0  can  be  expressed  as 


Py(P)  =  Hp  -  (C33  -  C44)S2<^’>  /i  Hp  -  (C33  -  C44)S^(P)' 


+  (Ci3  +  C44)S 


211/2 


P^(P)  =  (C12  +  C44)Sy<P)s2(P)  /J  Hp  -  (C33  -  C44)S 


+  (Ci3  +  C44)Sy(l’)S2^^^ 


2II/2 
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where  Hp  is  defined  as 

{2  2 
(Cii  -  C44)Sy(P)  +  (C33  -  C44)S2<P) 

rr  2  2I2 

+  j  (Cii  -  C44)Sy(P>  +  (C33  -  C44)S2(P) 

-  4Sy(P)\(P)^  (Cii  -  C44)(C33  -  C44) 

-  (Ci3  +  C44)2 

The  amplitude  ratios  of  the  reflected  SV  and  P  waves, 
and  ,  for  an  incident  SV  wave  of  unit  amplitude  whose  angle 

of  incidence  is  less  than  the  critical  angle  can  be  determined 
from  Eq.  (28).  For  a  given  incident  SV  wave,  the  values  of  the 
components  of  the  slowness  vector  along  the  y  and  z  axes,  Sy^^^ 
and  82^^^  and  the  values  of  the  components  of  the  unit  vector  of 
particle  displacement  along  the  y  and  z  axes,  Py^^^  and  P^^^^ 
are  defined  as  part  of  the  specification  of  the  incident  SV  wave. 
The  values  of  the  components  of  the  slowness  vector  of  the 
reflected  SV  wave  along  the  y  and  z  axes,  Sy^^^^  and  are 

determined  from  Eq.  (16)  and  (17),  respectively.  The  values  of 
the  components  of  the  unit  vector  of  particle  displacement  along 
the  y  and  z  axes  for  the  reflected  SV  wave,  Py^^^^  and  P2^^^\ 
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are  obtained  by  substituting  and  into  Eq.  (29). 

The  values  of  the  coefficients  for  in  Eq.  (28)  are  thus 

obtained.  Since  the  y-component  of  the  slowness  vector  of  the 
reflected  P  wave,  ,  is  equal  to  that  of  the  incident  SV 

wave,  Eq.  (16),  the  value  of  the  z-component  of  the  reflected  P 
wave,  is  obtained  by  substituting  Sy^^)  into  the  slowness 

surface  for  the  P  wave,  Eq.  (15).  On  substitution  of  the  values 
of  and  Sy^^)  into  Eq.  (30),  the  values  of  the  unit  vector 

of  particle  displacement  along  the  y  and  z  axes  for  the  reflected 
P  wave,  Py(P)  and  P2^^^,  are  determined.  The  values  of  the 
coefficients  for  in  Eq.  (28)  are  thus  obtained.  The  values 

of  the  amplitudes  of  the  reflected  SV  and  P  waves,  A^^^^  and 
are  then  obtained  by  solving  Eq.  (28)  with  the  thus- 
determined  values  of  the  coefficients  for  A^^^)  and  A^^)  in  Eq. 
(28)  for  a  given  incident  SV  wave.  By  varying  the  angle  of  inci¬ 
dence  of  the  incident  wave.  By  varying  the  angle  of  incidence  of 
the  incident  SV  wave  of  unit  amplitude  and  by  repeating  the 
procedures  described  above,  the  amplitude  ratios  of  the  reflected 
SV  and  P  waves  to  the  the  incident  SV  wave  are  obtained  as  func¬ 
tions  of  the  angle  of  incidence.  For  an  incident  SV  wave  in  the 
unidirectional  fiberglass  epoxy  composite  as  shown  in  Fig.  1,  the 
amplitude  ratios  of  the  reflected  SV  and  P  waves  to  the  incident 
SV  wave  versus  the  angle  of  incidence  are  plotted  in  Fig.  5  from 
zero  to  the  critical  angle  52®. 
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5 .  Existence  of  Reflected  Surface  Wave  for  Angle  of  Incidence 


Equal  to  or  Greater  Than  Critical  Angle 

It  is  seen  from  Fig.  5  that  at  the  critical  angle 
52°,  the  amplitude  ratio  of  the  reflected  SV  wave  to  the  incident 
SV  wave  is  minus  one  and  the  amplitude  ratio  of  the  reflected  P 
wave  to  the  incident  SV  wave  is  -2.98  instead  of  being  equal  to 
zero.  This  suggests  the  possibility  of  the  existence  of  a  sur¬ 
face  wave  travelling  parallel  to  the  plane  boundary  for  the  angle 
of  incidence  equal  to  or  greater  than  the  critical  angle.  In  the 
case  of  an  incident  SV  wave  in  a  semi- infinite  isotropic  medium 
with  an  angle  of  incidence  greater  than  the  critical  angle,  there 
exist  a  reflected  SV  wave  and  a  surface  P  wave  whose  amplitude 
decays  exponentially  with  the  perpendicular  distance  from  the 
plane  boundary  [4] . 

Consider  an  SV  wave  travelling  in  the  plane  x  =■  0  in  the 
semi -infinite  transversely  isotropic  medium  which  is  incident  on 
the  plane  boundary,  with  the  angle  of  incidence  equal  to  or 
greater  than  the  critical  angle.  Assume  that  there  exists  a  sur¬ 
face  wave  travelling  parallel  to  the  plane  boundary;  that  is,  the 
x-component  and  the  z -component  of  its  slowness  vector,  Sj^  and 
S^,  vanish,  whereas  the  y-component  of  its  slowness  vector,  Sy, 
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does  not  vanish,  as  shown  in  Fig.  6(a).  Let  the  reflected  sur¬ 
face  wave  be  represented  as 


(u,  V,  w)  =  A(^)exp(-C]^z)(Px^^\ 

exp{ia?(Sy^^)y  -  t) }  (31) 

where  Px^^^ »  components  of  a  unit  vector 

of  particle  displacement  of  the  reflected  surface  wave  along  the 
X,  y  and  z  axes,  respectively;  is  the  amplitude  of  particle 

displacement  of  the  reflected  surface  wave;  Sy^^^  is  the  y-com- 
ponent  of  the  slowness  vector  of  the  reflected  surface  wave;  and 
is  an  undetermined  real  constant  which  is  equal  to  or  greater 
than  zero. 

Now,  rearrange  Eq,  (31)  into  the  following  form: 

(u,  V,  w)  =  a(S)(Pj^(S)  _  Py(S)^ 

exp{iw[Sy(^V  +  i(Ci/w)z]  -t}  (32) 

On  substitution  of  Eq.  (32)  into  Eq.  (11),  we  find  that 

“  0 

-  C44(Ci/a))2  -  p]Py(S) 
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+  (Ci3  +  C44)Sy<S)i(c^/j^)P^(S)  .  0 


(Ci3  +  C44)Sy(S)i(C3^/^)P  (S) 


+  [C44S„(S)  .  C33(Ci/«)2  -  -  0 


The  condition  for  the  existence  of  non- trivial  solutions  for 
Py(S)  and  is  expressed  by  setting  the  determinant  of  the 

matrix  of  the  coefficients  of  ^nd  in  Eq.  (33)  equal 

to  zero: 


[CiiS„(S)  .  C44(Ci/a>)2  -  p]  (0^3  +  C44) Sy(S) i(Ci/w) 


(Ci3  +  G44)Sy(S)i(Ci/w)  [C44Sy<S)  -  C33(Ci/w)2  -  p] 


Expanding  Eq.  (34) ,  wt  obtain 


12  r  2 

*^33*^44  C33(p  -  ) 


2  2 

+  C44(p  -  C44Sy<S)  )  +  (Ci3  +  C44)2Sy(S)  (Ci/w)2 


2  2 

+  (p  -  CiiSy(S)  )(p  -  C44Sy(S)  )  .  0 


Eq.  (35)  is  a  quadratic  equation  in  (Cj^/w)^;  designate  its 
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coefficients  as  follows: 


A  =  C33C44 


2  2 
B  =  C33(p  -  CiiSy(S)  )  +  C44(p  -  C44Sy(S)  ) 

+  (Ci3  +  C44)2Sy(S)^  (36) 

2  2 

C  =  (p  -  CiiSy(S)  )(p  -  C44Sy(S)  ) 


It  is  apparent  from  Eq.  (36)  that  A  is  always  positive. 

By  definition,  at  the  critical  angle,  the  angle  of  reflection 
of  the  reflected  P  wave  is  equal  to  90®  [4] ;  accordingly, 

$2^^^  in  Eq.  (20)  is  equal  to  zero.  Therefore,  the  value  of  b  in 
Eq.  (16)  corresponding  to  the  critical  angle  is  thus  determined 
by  setting  and  in  Eq.  (15)  equal  to  zero,  and  is  found  to 
be  .  As  stated  following  Zq.  (19),  the  angle  of 

reflection  of  the  reflected  SV  wave  is  equal  to  90®  when  the 
SV  wave  is  tangentially  incident  on  the  plane  boundary;  accord¬ 
ingly,  in  Eq.  (18)  is  equal  to  zero.  Therefore,  the  value 

of  b  in  Eq.  (16)  corresponding  to  the  incident  SV  wave  whose 
slowness  vector  is  parallel  to  the  plane  boundary  and  is  in  the 
plane  x  =  0  is  determined  by  setting  and  S2  in  Eq.  (14)  equal 
to  zero,  and  is  found  to  be  ,  The  necessary  condition 

for  the  occurrence  of  the  critical  angle  phenomenon  for  an  SV 
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wave  travelling  in  the  plane  x  =  0  is  that  the  value  of 
(p/C44)^/2  has  to  be  greater  than  that  of  [5]  Xake 

numerical  values  of  and  (p/Cn)^/^  for  the 

unidirectional  fiberglass  epoxy  composite  shown  in  Fig.  1  from 
[1]  as  follows:  =  0.418 (s/km)  and  (^7044)^/2  = 

0.647 (s/km),  as  shown  in  Fig.  3.  Therefore,  these  values  satis¬ 
fying  the  necessary  condition  on  the  occurrence  of  the  critical 
angle  phenomenon  for  an  incident  SV  wave  in  the  unidirectional 
fiberglass  epoxy  composite. 

It  has  been  shown  [5]  that  when  the  angle  of  incidence  of  an 
incident  SV  wave  travelling  in  the  plane  x  =  0  is  equal  to  or 
greater  than  the  critical  angle,  the  value  of  b  in  Eq.  (16)  lies 
between  (p/C]^;^)^^^  (p/C44)^/^,  that  is;  (p/C^^l)^/^  <  b  < 

(p/C44)V2^  as  shown  in  Fig.  3.  In  addition,  it  follows  from  Eq. 
(5)  that  the  y- component  of  the  slowness  vector  of  the  reflected 
surface  wave,  is  equal  to  b  in  Eq.  (16).  As  a  result,  the 

value  C  in  Eq.  (36)  is  equal  to  or  less  than  zero  when  the  angle 
of  incidence  of  an  incident  SV  wave  travelling  in  the  plane  x  -  0 
is  equal  to  or  greater  than  the  critical  angle ;  that  is ,  C  <  0 . 
Accordingly,  the  value  of  -  4AC  defined  in  Eq.  (36)  is  equal 
to  or  greater  than  zero  for  the  angle  of  incidence  equal  to  or 
greater  than  the  critical  angle;  that  is. 
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b2  -  4AC  =  b2  +  4a|c|  >  0 


(37) 


Thus,  the  possible  solutions  to  the  quadratic  equation  in 


in  Eq .  (35) .are 


,  -B  ±  (b2  -  4AC)V2 

(Ci/«)2 - ^ - i - 


However,  it  follows  from  Eq.  (37)  that  the  value  of  B^  -4AC  is 
equal  to  or  greater  than  that  of  B^  for  any  angle  of  incidence 
equal  to  or  greater  than  the  critical  angle.  Consequently,  the 
solution  to  (Ci/w)2,  from  Eq.  (38),  is 


(Ci/c)^  = 


2  _  -B  +  (b2  -  4AC)  V2 


Accordingly,  there  exists  a  real  constant  Ci  given  in  Eq.  (31) 
when  the  angle  of  incidence  of  the  incident  SV  wave  travelling  in 
the  plane  x  =  0  is  equal  to  or  greater  than  the  critical  angle, 
and  is  expressed  from  Eq.  (39),  as 


-B  +  (b2  -  4AC)V2  ■ 

Cl  =  w  - - -  >  0 

2A 


The  existence  of  a  real  constant  Cj^  in  Eq.  (40)  confirms  the 
existence  of  the  surface  wave  assumed  in  Eq.  (31). 


It  is  therefore  concluded  that  when  the  angle  of  incidence  of 
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an  incident  SV  wave  travelling  in  the  unidirectional  fiberglass 
epoxy  composite  as  shown  in  Fig.  1  is  equal  to  or  greater  than 
the  critical  angle,  only  an  SV  wave  is  reflected  back  into  the 
medium,  and  the  reflected  P  wave  degenerates  into  a  surface  wave 
which  travels  parallel  to  the  plane  boundary  and  whose  amplitude 
decreases  exponentially  with  perpendicular  distance  form  the 
plane  boundary,  as  shown  in  Fig.  6(b). 

6 .  Balance  in  Energy  Flux  Normal  to  Plane  Boundary 

The  balance  in  energy  flux  normal  to  the  plane  boundary  must 
be  satisfied  [2] .  For  the  angle  of  incidence  less  than  the  crit¬ 
ical  angle,  the  balance  in  energy  flux  normal  to  the  plane  bound¬ 
ary  z  =  0,  as  shown  in  Fig.  2,  is  expressed  as  [2] 

F^d)  +  F2(SV)  =  0  (41) 

where  F2d),  z-components  of  the  energy 

fluxes  of  the  incident  SV  wave,  the  reflected  SV  wave  and  the 
reflected  P  wave,  respectively. 

The  z -component  of  the  energy  flux  of  an  incident  SV  wave  of 
unit  amplitude  travelling  in  the  plane  x  -  0  is  [6] 
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(42) 

+  C44Py(I)p,(I)Sy(I)  +  C33P,(I)'s,(I)) 

Similarly,  the  z-components  of  the  energy  fluxes  of  the  reflected 
SV  and  the  reflected  P  waves  are  [6] 

F^<r)  -A(I')V<C44Py(P)^S^(f)  f  C13P  (P)p^(P)Sy(P) 

(43) 

+  C44Py(P>P,(P)Sy(P>  *  C33P,<P)^S,<P)) 

and 

F^(SV)  .A(SV)2„2(c44Py(SV)2s^(SV)  Ci3Py(SV)p^ (SV)s^(SV) 

+  C44Py(SV)p^(SV)s^(SV)  H.  C33P^<SV)2s^(SV)) 

Since  the  frequency  term  is  common  to  Eqs.  (42),  (43)  and 
(44) ,  the  balance  in  energy  flux  normal  to  the  plane  boundary, 

Eq.  (41),  is  not  affected  by  assuming  the  value  of  the  frequency 
(j>  to  be  equal  to  unity.  Accordingly,  subsequent  calculations  of 
the  values  of  the  z-components  of  the  energy  fluxes  of  the  inci¬ 
dent  SV  wave  and  the  reflected  SV  and  P  waves,  and 

Fz^^K  ate  done  by  assuming  the  radian  frequency  u  to  be  equal  to 


one . 


The  value  of  the  z-components  of  the  energy  fluxes  in  Eqs . 

(42)  ,(43)  and  (44)  are  obtained  similarly  to  the  calculations  of 

the  amplitude  ratios  of  the  reflected  P  and  SV  waves  to  the  inci¬ 
dent  SV  wave  for  the  angle  of  incidence  less  than  the  critical 
angle.  For  a  given  incident  SV  wave  of  unit  amplitude,  the 
y-components  of  the  slowness  vector  and  of  unit  vectors  of  par¬ 
ticle  displacement,  Py^^^^  and 

Py(P) ,  and  the  z- component  of  the  slowness  vectors  and  of  the 
unit  vectors  of  particle  displacement,  82^^^  S2^^\ 

P  and  P_^^)  of  the  incident  SV  wave  and  the  reflected 

SV  and  P  waves  are  determined  first.  Then  combining  the  ampli¬ 
tudes  of  the  reflected  SV  and  P  waves,  and  with  the 


z-components  of  the  energy  fluxes  of  the  incident  SV  wave 


and  the  reflected  SV  and  P  waves, 

thus  obtained  from  Eqs.  (42),  (43)  and  (44). 

For  an  SV  wave  of  unit  amplitude  travelling  in  the  plane  x  = 

0  in  the  unidirectional  fiberglass  epoxy  composite  shown  in  Fig. 

1  incident  on  a  plane  boundary,  the  z-components  of  the  energy 
fluxes  of  the  reflected  SV  and  P  waves  and  the  energy  flux  of  the 
incident  SV  wave  are  plotted  in  Fig.  7  from  zero  to  the  critical 
angle  52°,  with  the  value  of  the  frequency  w  in  Eqs.  (42),  (43) 


33 


and  (44)  equal  to  one.  The  balance  in  energy  flux  normal  to  the 
plane  boundary,  Eq.  (41),  is  also  checked  and  is  shown  in  Fig.  7. 

When  the  angle  of  incidence  of  an  incident  SV  is  equal  to  or 
greater  than  the  critical  angle,  only  an  SV  wave  is  reflected 
back  into  the  medium  with  its  amplitude  ratio  to  the  incident  SV 
equal  to  minus  one,  and  the  reflected  P  wave  degenerates  into  a 
surface  wave  whose  energy  flux  is  parallel  to  the  plane  boundary 
[6].  As  a  result,  the  energy  flux  normal  to  the  plane  boundary 
due  to  the  surface  wave  is  equal  to  zero.  The  balance  in  energy 
flux  normal  to  the  plane  boundary  z  —  0,  as  shown  in  Fig.  6(a), 
for  the  angle  of  incidence  equal  to  or  greater  than  the  critical 
angle  is  expressed  as  [2] 

f(I)  +  =  0  (45) 


The  relationships  between  the  components  of  the  unit  vector 
of  particle  displacement  of  an  incident  SV  wave  travelling  in  the 
plane  x  -  0  along  the  y  and  z  axes,  those  of 

the  reflected  SV  wave  can  be  obtained  by  substituting  Eqs.  (16) 
and  (17)  into  Eq.  (29)  and  can  be  expressed  as 


(46) 
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On  substitution  of  Eqs.  (16),  (17)  and  (46)  into  Eqs.  (42) 


and  (44) ,  it  is  found  that 

F^CI)  .  .F^<SV) 


(47) 


Accordingly,  the  balance  in  energy  flux  normal  to  the  plane  bou¬ 
ndary,  Eq  (45),  is  satisfied  when  the  angle  of  incidence  of  an 
incident  SV  wave  is  equal  to  or  greater  than  the  critical  angle. 
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ACOUSTO-ULTRASONIC  NON^CONTACT  INPUT>OUTPUT  CHARACTERIZATION  OF 
UNIDIRECTIONAL  FIBERGLASS  EPOXY  COMPOSITE  PLATE 

It  has  been  shown  that  the  unidirectional  fiberglass  epoxy 
composite  shown  in  Fig.  1  may  be  modelled  as  a  homogeneous  trans¬ 
versely  isotropic  continuum.  For  the  axes  shown  in  Fig.  1,  the 
isotropic  plane  of  the  equivalent  continuum  lies  in  the  midplane 
of  the  plate  [1].  A  cartesian  coordinate  sysfem  (x,  y,  z)  is 
chosen  so  that  the  x-y  plane  is  the  isotropic  plane;  thus,  the 
upper  and  lower  faces  are  at  z  «  h/2  and  z  *=  -h/2,  respectively, 
where  h  is  the  plate  thickness.  The  properties  of  the  equivalent 
continuum  model  of  the  unidirectional  fiberglass  epoxy  composite 
plate  to  be  considered  are  [1] 

h  =  0 . 1  m 

Cii  =  10.581  X  10^  N/m2 

Ci3  =  4.67  X  10^  N/m^  (48) 

C33  -  40.741  X  10^  N/m2 
C44  =  4.422  X  10^  N/m2 
Cgg  =  3.243  X  10^  N/m2 
p  =  1850  kg/m^ 

Non- contact  transmitting  and  receiving  transducers  are 
located  on  the  same  face  of  a  fiberglass  epoxy  composite  plate 
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specimen,  as  shown  in  Fig.  8.  The  unidirectional  fiberglass 
epoxy  composite  plate  specimen  is  considered  as  a  plate  of  thick¬ 
ness  h  and  of  infinite  planar  (x-y)  extent.  The  input  electrical 
voltage  to  the  transmitting  transducer  is  V^(t)  and  the  output 
electrical  voltage  from  the  receiving  transducer  is  V^Ct)  where  t 
represents  time.  The  transmitting  transducer  converts  an  input 
electrical  voltage  into  a  stress,  whereas  the  receiving  trans¬ 
ducer  converts  a  displacement  associated  with  stress  waves  trav¬ 
elling  in  the  plate  into  an  output  electrical  voltage.  In  the 
following  analysis,  only  the  SV  waves  are  traced.  The  SV  waves 
which  are  generated  by  the  transmitting  transducer  located  above 
point  0  experience  multiple  reflection  at  each  face  of  the  plate, 
and  then  reach  the  receiving  transducer  located  above  point  M,  as 
shown  in  Fig.  9.  Since  the  isotropic  plane  lies  in  the  midplane 
and  is  parallel  to  both  the  top  and  bottom  faces  where  the  mul¬ 
tiple  reflections  occur,  the  angle  of  reflection  of  the  reflected 
SV  wave  is  equal  to  the  angle  of  incidence  of  an  incident  SV  wave 
for  each  reflection  at  each  face  of  the  plate.  Accordingly,  the 
SV  waves  travelling  from  the  input  0  to  the  output  M  may  be  con¬ 
sidered  as  waves  propagating  in  a  semi- infinite  transversely  iso¬ 
tropic  medium  and  travelling  to  point  M'  as  if  there  were  no 
bottom  face,  as  shown  in  Fig.  9. 
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1 •  Delay  Time  and  Phase  Velocity 


Let  the  input  0  and  the  output  M  lie  in  the  y-z  plane. 

Assume  the  numbet  of  reflections  at  the  bottom  face  experienced  by 
the  SV  waye  in  trayelling  from  the  input  0  to  the  output  M  is  n, 
as  shown  in  Fig.  9.  With  respect  to  the  z  axis,  the  angle  of 
incidence  of  the  SV  waye  at  each  face  of  the  plate  is  8 ,  and  the 
total  distance  travelled  by  the  wave  is  R^.  From  the  geometry  in 
Fig.  9, 


8  =  tan'^(i/2nh)  (49) 

where  i  is  the  separation  distance  between  the  input  0  and  the 
output  M, 

Rn  =  £/sin8  (5(» 

The  delay  time  tj^  for  the  SV  wave  to  reach  the  receiving  trans¬ 
ducer  is 


t^n  =  Rn/Ci(^)  (51) 

where  C^CS)  is  the  directionally- dependent  phase  velocity  of  the 
SV  wave.  The  phase  velocity  of  an  SV  wave  in  the  unidirec¬ 
tional  fiberglass  epoxy  composite  plate  is  [1] 

Ci(.8)  =  [(C44  +  C^i^sin^^  +  033003^^  -  yE)/2/p]  (52) 
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where 


E  =  [(Cii  -  C44)sin2^  +  (C44  -  033)003^^] 

+  4(Cj^3  +  C44)^sin^^cos^S  ; 

^11  >  ^13  >  ^33'  *^44  P  given  by  Eq.  (48). 

The  phase  velocity  as  a  function  of  the  angle  of  incidence  6 
is  shown  in  Fig.  10. 

The  delay  time  is  then  computed  when  the  number  of  reflec¬ 
tions  n  at  the  bottom  face  of  the  plate  is  equal  to  10,  100,  300 
or  500.  The  numerical  results  are  shown  in  Fig.  11  where  the 
delay  time  tjj  is  plotted  as  the  ordinate,  and  the  dimensionless 
separation  i/h  is  plotted  as  the  abscissa,  for  values  of  zero  to 
300. 

2 .  Displacements  Detected  bv  Receiving  Transducer 

The  displacements  detected  by  the  non-contact  receiving 
transducer  above  point  M,  radiated  by  the  non-contact  transmit¬ 
ting  transducer,  are  assumed  to  be  equivalent  to  the  displacement 
at  point  M'  associated  with  the  SV  wave  propagating  in  a  semi¬ 
infinite  transversely  isotropic  medium.  Except  for  the  reflection 
coefficients  at  each  face  (to  be  discussed  later) ,  the  displace¬ 
ment  is  computed  as  if  there  were  no  bottom  boundary,  as  shown  in 
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Fig.  9.  The  displacement  at  point  M'  is  approximated  by  the 
far- field  asymptotic  solution  for  large  of  an  infinite  trans¬ 
versely  isotropic  medium  subjected  to  a  harmonic  point  force. 

Consider  an  infinite  transversely  isotropic  medium  in  which 
the  z  axis  of  a  rectangular  cartesian  coordinate  system  0(x,  y, 
z)  is  the  zonal  axis  of  the  medium,  and  the  x-y  plane  coincides 
with  isotropic  plane,  as  shown  in  Fig.  12.  The  equations  of 
motion  including  the  body  force  are  [7] 

^xx,x  ^xy,y  ^xz,z  +  pX  =  pu,tt  (53) 

^xy,x  ■*"  ^yy,y  ^yz,z  +  pY  =  pv,tt  (54) 

^xz,x  ^yz,y  ^zz,z  +  pZ  pw,tt  (55) 

where  (r,  s  =  x,  y  and  z)  are  the  normal  (r  =  s)  and  shear 

(r  5^  s)  stresses  with  respect  to  the  chosen  coordinate  system 
0(x,  y,  z) ;  u,  V  and  w  are  displacement  components  of  a  point  in 
the  medium  along  the  x,  y  and  z  axes,  respectively;  X,  Y  and  Z 
are  the  components  of  the  body  force  along  the  x,  y  and  z  axes, 
respectively;  p  is  the  density;  t  is  time  ;  and  denotes 
partial  differentiation  with  respect  to  the  variable  which  fol¬ 
lows  . 
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Combining  Eqs.  (1)  and  (55)  gives 


r  tt  =  (^13  ^44_^  ^  (-^^)  r.xx  +  (  )  r.yy 

p  P  P 

C33 

+  ( — —)  r,zz  +  z,z 

p 


where  T  and  A  are  given  by  [6] 


r  =  w,z  and  A  =  u,x  +  v,y 


By  differentiating  Eq.  (53)  with  respect  to  x  and  Eq.  (54)  with 
respect  to  y,  we  find,  on  addition  of  the  resulting  equations  and 
using  the  appropriate  stress-strain  relation,  Eq.  (1),  that 


Cl  'x  +  Caa  ^44 

A,tt  - —)  (r,xx  +  r,yy)  + -  A,zz 

P  P 


+  (A, XX  +  A,yy)  +  X,x  +  Y,y 

P 


For  a  harmonic  point  force  at  the  origin,  the  body  forces  may 
be  taken  of  the  form  [7] 


X  =  Xo5(x)S(y)5(z)e- 


Y  =  Yo5(x)S(y)S(z)e- 


Z  -  ZQS(x)5(y)6(z)e“ 


where  5(r)  (r  =  x,  y  and  z)  is  the  Dirac  delta  function. 


Express  F  and  A  as  threefold  Fourier  Integrals  [7]  as 


follows : 


r(x,y,z,t)  =  r(Sx,Sy.,S2,t) 


exp{iw(Sj^x  +  SyY  +  S^z  -  t)  jdS^^dSydS^ 


A(x,y,z,t)  «  A(S5^,Sy,S2,  t) 


exp{iw(S^x  +  SyY  +  S^z  -  t))dS^dSydS2 


where 


r(Sx,Sy,S2,  t)  -  l/Sn'^  r(x,y,z,t) 


exp{  icoC-Sj^x  -  Syy  -  S^z  +  t)}dxdydz. 


A(Sx,Sy,S2,t)  =  I/Sjt^  A(x,y,z,t) 


exp{iw(-SxX  -  Syy  -  S^z  +  t)}dxdydz. 
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Similarly,  X,x,  Y,y  and  Z,z  can  be  expressed  as  threefold  Fourier 
integrals 


X,x  =  iwS^  X  exp { icj(Sj^x  +  Syy  +  ■  t))dS^dSydS2 


Y,y  —  icoSy  Y  exp{iw(S2^x  +  Syy  +  ^2^  "  jdS^dSydS^ 


Z,z  =  10)82  ^  exp{io)(SxX  +  Syy  +  S2Z  -  t))dS2dSydS2 


where 


X  =  l/8jr^  Xo5(x)5(y)5(z)e- 


exp{io;(-SjjX  -  Syy  -  S2Z  +  t)dSxdSydS2  =■  Xq/Sw 


and,  similarly.  Y  =  Yq/Stt^;  Z  -  Zo/8jr-^. 
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In  fact,  H(Sjj,  Sy,  Sg)  =  0  represents  the  two  sheets  of  slowness 
surface,  one  for  a  P  wave  and  one  for  an  SV  wave  [8].  As  a 
result  of  Eqs.  (63)  and  (64), 
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r(x,y,z,t)  = 


(Ci3  +  C44)/p 


85r^c<)H(S^,  Sy ,  S2) 


(iSjjXo  +  iSyYg) 


fCll/p(Sx^  +  Sy2)  +  033/^82^  -1' 
87r^a>H(Sjj,  Sy ,  S2)  , 

exp{i«(SxX  +  SyY  +  S2Z  -  t) ) 


iSzZo 


(66) 


A(x,y,z,t)  = 


C44//>(Sx^  +  +  C33/pS2^  -1 


8w^H(Sjj,  Sy ,  S2) 

2  _i.  Q  2y 


(iSxXo 


f-(Ci3  -f  C44)/p  (Sx  +  Syjl  . 
[  85r^coH{Sjj,  Sy,  S2)  J 


(67) 


exp{iw(SxX  +  Syy  +  S2Z  -  t) } 


The  asymptotic  solution  at  a  large  distance  from  the  point 
force  is  obtained  by  applying  the  theory  of  residues,  the  method 
of  stationary  phase,  and  the  radiation  condition  [8]  as 


r(x,y,z,t) 


w  N 

~  ;r^  2-^^n 

2jrR  n=l 


-(C13+C44)  ^  2,,^  V  V  ^ 


.  1^11  2  .  s  2.  C33  2 

+  I -  (Sxn  +  Syn  7  +  ^zn 

P  P 


[ 

1^“^  '“^xn  ^  “yn 

exp{iw(SxnX  +  S^^y  +  82^2  -  t) } 


^^zn^o 


(68) 
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CO  N 

A(x,y,z,t) - S  A_A„  i 

2^rR  n=l  ' 


—  (Sxn^+  S  Szn^-  1 

/I  ''  r\ 


(ISj^n^o'^^SynYo)  + 


-(C13+C44) 


^yn  ) 


} 


exp{iw(SxnX  +  Sy^y  +  S^nZ  -  t) ) 


(69) 


where  R  is  the  distance  from  the  origin  0  where  the  point  forces 
are  applied  to  the  location  of  interest  Q  in  the  medium,  as  shown 
in  Fig.  12;  Sy^,  are  points  on  the  slowness  surface 

for  an  SV  wave  where  the  normal  is  parallel  to  the  OQ  direction; 

N  is  the  total  number  of  points  Sy^^,  S^n)  are  points  over 

which  the  summation  must  be  performed  ;  is  the  amplitude 
coefficient  and  is  given  by 


A 


n 


+  H,Sy2  + 

Kl 


1/2 


(70) 


where  is  given  as 


kjj  =  S  [H,S22(H,SxSx  H.SySy  -  H,SxSy2) 

+  2H,Sx  H,Sy(H,SxS2  H,SyS2  -  H.S^Sy  H.S^S^); 


(71) 


where  S  denotes  the  sum  with  respect  to  cyclic  permutation  of  Sj^, 
Sy  and  S2;  the  symbol  |  |  denotes  "the  magnitude  of";  A^^  is 
evaluated  at  points  (S^n.  Sy„,  S2n)  on  the  slowness  surface  for 
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an  SV  wave  where  the  normal  is  parallel  to  the  OQ  direction, 
exclvisive  of  those  singular  points  ^yn*  ^zn^  resnlting  in 

in  Eq.  (71)  equal  to  zero;  and  is  the  phase  constant  and  is 
determined  as  follows:  =  1  if  >  0  or  A^  =  i  if  <  0. 


The  displacement  components  along  the  x,  y  and  z  axes,  u,  v 
and  w  due  to  an  SV  wave  can  be  obtained  by  direct  integration  of 


the  definition  of  T  and  A  in  Eq.  (56),  and  are  given  as  [8] 


A 

- R - R  A 

+  Sy^) 


V  =  - R^- - ?5  A  (72) 

iw(Sx^  +  Sy'^) 

1 

w  =  -  r 

Substitution  of  Eqs .  (68)  and  (69)  into  Eq.  (72)  gives  the 
asymptotic  solutions  at  a  large  distance  to  the  displacement  com¬ 
ponents  along  the  x,  y  and  z  axes,  u,  v  and  w,  as  follows: 


VAX  f 

2;rR  n=l^  "  t  (^xn^  ^ 


S 

xn  ’  ^yn  ' 


/c  2.C  S  2  _  ^ 

IP  P 


"zn 


(iSxn^Xo+iSyn^Yo)  +  +  Sy^^)  iS^^^ 


exp{it<>(SxnX  +  S^^y  +  -  t)  }  (73) 
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-i  N  S 

2^R  n=l^  ^ 


— ^  ( g  ^  +  S  I  ^ ^  ^  1 

^^xn  ^  ^yn  '  ^  ^zn  “  ^ 

P 


(lS„2x„+ISy„2Y„)  +  — (S^„2  +  S  2)  iS^„Z„ 


exp{iw(SxnX  +  Sy^y  +  -  t) 


1  N 

2jrR  n=l  ^  ^ 


-(C13+C44) 


^zn^  ^SjjjjXo+iSynYo) 


exp{ia)(SxnX  +  Sj^y  +  S^nZ  -  t) } 


3 .  Directivity  Functions 


The  shear  stress  ty^  and  the  normal  stress  associated 
with  the  SV  waves  reaching  the  point  M'  in  Fig.  9  are  used  to 
study  their  associated  directivity  functions.  The  asymptotic 
shear  stress  Ty^  and  the  asymptotic  normal  stress  are  given 
^y  substituting  Eqs.  (68)  through  (72)  into  Eq.  (1)  as 


N  f  [■  (^13'*’^44)  0 

^JiVn  |C33|^ - - -  S^„2(is^^X„+iSy„Y„! 


Cii  ^ 

+  - (S™^ 


s  2 

^zn 

P 


^vn  )  ”*■  ^zn  ‘  ^K^zn^o 
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^13 


^44, „  2  a.  c  2^  .  ^33  2 


LI  P 


■^^xn  ^  ^yn  ^  ■*" 


L)(i 


^zn  “  ^ 


*  (S,„2  +  Sy,2)  lS,„zJ. 

6xp { lco( +  ^yn.y  ^ 


(76) 


~  2xR 


jsyn  * 


-(C13+C44) 


'zn 


( ^  ^  xn^o"*"  ^  ^  yn^o  ^ 


+  (— (S^n^  +  +  — 


Szn^  •  DlZo 


zn 


c  2  c  2 

^xn  ^  ^yn 


[ -  (^xn^  ^  ^zn^  "  1] 


+  (Sxn^  +  S  2)  is^^Zo 

P 


exp { iw( Sjjj^x  +  Syj^y  +  -  t) ) 


(77) 


The  directivity  functions  associated  with  the  normal  stress 
.r2z  in  Eq.  (76)  and  the  shear  stress  Ty^  in  Eq.  (77)  will  be 
evaluated  for  the  case  of  =  0.  Due  to  the  axial  symmetry  with 
respect  to  the  zonal  axis,  the  z  axis,  of  the  transversely  iso¬ 
tropic  medium,  the  values  of  the  directivity  functions  thus 
obtained  hold  for  all  values  of  Sjj. 


Consider  the  case  of  the  point  force  acting  along  the  y 


49 


direction  only;  that  is,  Yq  0  but  Xq  =  Zq  =■  0  in  Eq.  (58). 

Y 

The  directivity  function  Dy^  of  the  shear  stress 
associated  with  an  SV  wave  is  evaluated  from  Eq.  (77)  along  the 
arc  y^  +  =  1  in  the  positive  y-z  quadrant  by  setting  Yq  =  R  = 

1; 


2k 


N  "*^13  9  C33  „ 

^^n^zn^  ^yn  "*■  ^zn 

n=l  P  p 

exp{ia>(Sy^y  +  S^n^)) 


1) 


(78) 


where  the  summation  S  is  done  vectorially  in  a  complex  plane  over 
each  type  of  wave  passing  a  given  point  (0,  y,  z)  on  the  arc  y^  + 

O 

z^  =  1;  each  wave  corresponds  to  a  point  (0,  Sy,^,  S^^)  on  the 
slowness  surface  of  an  SV  wave  at  which  the  normal  is  parallel  to 
the  radius  vector  of  the  given  point  (0,  y,  z) .  Similarly,  the 

Y 

directivity  function  D^z  of  the  normal  stress  Tzz  is,  from  Eq. 
(76),  for  Yq  =  R  -  1; 


N 

S 

n=l 


^^n^yn 


C13C44  „  ? 
7“  V 


C33C44  „  0 

^zn 


C 


13 


exp{iw(Syj^y  +  S^nZ)) 


(79) 


Next,  consider  the  case  of  the  point  force  acting  along  the  z 
direction  only;  that  is,  Zq  f*  0  but  Xq  -  Yq  -  0  in  Eq.  (58). 

Z 
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The  directivity  function  Dy^  of  the  shear  stress  Ty^ 
associated  with  an  SV  wave  is  evaluated  from  Eq.  (77)  along  the 
arc  y^  +  =  1  in  the  positive  y-z  quadrant  by  setting  Zq  =*  R  = 

1; 


^^55 

Ztt 


N 

V'^n^yn 
n=l  ^ 


+  ^33  “  ^44  -  ^13 
P 


exp{iw(Sy^y  +  S^nZ) ) 


(80) 


Similarly,  the  directivity  function  D22  of  the  normal  stress  r 
is,  from  Eq.  (76),  for  Zq  =  R  1; 


'zz 


CO 


Ztt 


N 

W^zn 

n-=l 


C11C33  -C13C44  -Ci3  ^  2 


■  ^33  .  2  r 

^zn  ^33 
P 


exp{iw(Syny  +  S^nZ)) 


(81) 


Finally,  consider  the  case  of  the  point  force  acting  along  the 
X  direction  only;  that  is,  Xq  0  but  Yq  =  Zq  =  0  in  Eq.  (58).  It 
follows  from  Eqs.  (76)  and  (77)  that  the  values  of  the  shear 
stress  ry2  and  the  normal  stress  associated  with  an  SV  wave 

travelling  in  the  plane  x  *  0  are  equal  to  zero. 


The  polar  diagrams  for  the  directivity  functions  of  the  shear 


stress 


and  the  normal  stress 


zz 


(D. 


yz  » 


-'ZZ  » 


Dy2  and  D^z  given 
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in  Eqs .  (78)  through  (81)  )  associated  with  the  SV  waves  propa¬ 
gating  in  the  unidirectional  fiberglass  epoxy  composite  as  shown 
in  Fig.  1  are  obtained  by  substituting  Eq.  (48)  into  Eqs.  (78) 
through  (81)  at  frequencies  of  0.75,  1.50  and  2.25  MHz. 

Numerical  results  are  shown  in  Figs.  13  through  24  where  the 
angle  of  incidence  9  given  in  Eq.  (49)  is  used  to  determine  the 
direction  in  which  the  values  of  points  (0,  Sy^,  S^n)  on  the 
slowness  surface,  where  the  normal  is  parallel  to  the  given 
direction,  is  thus  obtained. 

4 .  Assumptions  of  Transducers 

The  non-contact  transmitting  transducer  in  Fig.  8  is  assumed 
to  transform  an  electrical  voltage  into  a  uniform  stress;  how¬ 
ever,  the  non- contact  receiving  transducer  in  Fig.  8  transforms  a 
displacement  into  an  electrical  voltage.  The  approach  below  is 
similar  to  that  given  in  [9].  Referring  to  Fig.  8,  if  an  input 
voltage  of  amplitude  of  V  and  frequency  w  is  applied  according  to 

Vi(t)  =  Ve'^“^  .  (82) 

The  stress  a  that  is  introduced  into  the  specimen  plate  by  the 
non- contact  transmitting  transducer  is 


52 


<7(t)  - 


(83) 


where  Fj^Cw)  is  the  transduction  ratio  for  the  non-contact 
transmitting  transducer  in  transforming  a  voltage  to  a  stress 
and  (f>i  is  a  phase  angle.  In  Eqs.  (82)  and  (83),  the  complex 
harmonic  character  of  the  signals  is  expressed  in  the  complex 
notation  where  i  —  -y-l  and  only  the  real  part  of  these  and 
subsequent  equations  should  be  considered.  Thus,  the  amplitude  T 
of  the  applied  load  is  defined  as 

T  =  F]^(w)V  (8^) 

Similarly,  if  a  stress  wave  producing  a  displacement  d  of 
amplitude  D  and  frequency  w  that  is  detected  by  the  non-contact 
receiving  transducer  is  defined  as 

d(t)  -  (85) 

the  output  voltage  from  the  non-contact  receiving  transducer,  see 
Fig.  8,  is 

VQ(t)  =  F2(w)De‘^(“^  ‘^2^  (86) 

where  F2(w)  is  the  transduction  ratio  for  the  non- contact 
receiving  transducer  in  transforming  a  displacement  to  a  voltage, 
and  <t>2  is  a  phase  angle.  Thus,  the  amplitude  V'  of  the  output 
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electrical  voltage  is 


V'  =  F2(«)D  , 
The  characteristics  of  F]^(w)F2(w)  are  unknown  except  that  the 
dimensions  of  the  product  Fj^(w)  are  [kg/m^»sec^] . 


5 •  Steady-State  Output  Voltage  Amplitude  due  to  Multiple  Wave 
Reflections  In  Plate 


Since  the  SV  wave  traced  in  the  unidirectional  fiberglass 
epoxy  composite  plate  specimen  shown  in  Fig.  9  is  travelling  in 
the  y-z  plane,  it  follows  from  Eqs.  (73),  (74)  and  (75)  that  only 
the  displacement  components  along  the  y  and  z  axes,  v  and  w,  are 
detectable  at  the  point  M' .  Consider  first  the  point  force 
acting  along  the  y  direction  only;  that  is  Yq  0  but  Xq  =  Zq  =  0 
in  Eq.  (58). 

The  amplitude  of  the  y- component  displacement  evaluated  at 
the  point  M'  can  be  obtained  from  Eqs.  (50)  and  (74)  as 


nY  _  ^l(^yn>  ^zn)  „ 

^v  ‘T -  Iq 

^n 


where 


-l^^yn’  ®zn^  1/2k  +  (033/^)82 


exp{iw(Syny  +  S^nZ) 
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and  is  evaluated  along  the  arc  -  1.  Similarly,  the 

Y 

amplitude  of  the  z- component  displacement  evaluated  at  the 
point  M'  can  be  obtained  from  Eqs.  (50)  and  (75)  as 


where 


^2^^vn»  ^zn^  “  l/27r  S  ^44 Vp] SynS2:n 

n— 1 


exp(iw(Syj^y  +  S2nZ)) 

Next,  consider  the  point  force  acting  along  the  z  direction 
only:  that  is  Zq  ^  0  but  Xq  -  -  0  in  Eq.  (58). 

z 

The  amplitude  of  the  y- component  displacement  evaluated  at 
the  point  M'  can  be  obtained  from  Eqs.  (50)  and  (74)  as 


pZ  „  ^2(V»  Szn)  2^  (90) 

Similarly,  the  amplitude  of  the  z-component  displacement 
evaluated  at  the  point  M'  can  be  obtained  from  Eqs.  (50)  and  (75) 


.Z  ^3 


^3(^vn>  ^zn^ 


where 


f3(Syn’  =  V2jr  (Cii/p)Syn^  +  (033/^)82^^  -  1] 
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exp{iw(Syny  +  S^nZ)) 


Finally,  consider  the  point  force  acting  along  the  x  direc¬ 
tion  only;  that  is,  Xq  0  but  Yq  =  Zq  =  0  in  Eq.  (58).  It 
follows  from  Eqs.  (74)  and  (75)  that  the  amplitude  of  the  y-com- 
ponent  displacement  and  the  amplitude  of  the  z -component  dis¬ 
placement  vanish  at  the  point  M'  .  This  is  due  to  the  fact  that 
the  SV  wave  is  travelling  in  the  y-z  plane. 


According  to  Eqs.  (88)  through  (91),  the  amplitude  of  a  dis¬ 
placement  component  at  the  point  M'  ,  denoted  as  ,  can  be 
expressed  in  the  following  form  as 


= 


fi(S 


yn* 


Rn 


(92) 


where  T  is  the  applied  point  force  and  is  equivalent  to  either 
or  Zq  in  Eqs.  (88)  through  (91);  and  fi(Sy^,  S^n)  (i  =*  1,  2  or  3) 
is  determined  from  one  of  Eqs.  (88)  through  (91),  depending  on 
which  displacement  component  is  measured  and  along  which  direc¬ 
tion  the  point  load  is  applied. 


However,  with  the  bottom  boundary  present,  the  wave  is 
reflected  a  total  of  (2n  -1)  times,  as  shown  in  Fig.  9.  Thus, 
the  amplitude  of  displacement  at  the  point  M  is  and  is 
expressed  as 
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(93) 


%  =  Qsvsv  % 


where  Qgvsv  amplitude  ratio  of  the  reflected  SV  wave  to 

the  incident  SV  wave,  as  shown  in  Fig.  5,  and  is  a  function  of 
the  angle  of  incidence  0  defined  in  Eq.  (49).  Here,  we  have 
ignored  the  effects  of  mode  conversion. 


The  amplitude  of  the  output  voltage  form  the  non-contact 
receiving  transducer  is  V'  and  can  be  obtained  by  substituting 
Eqs.  (92)  and  (93)  into  Eq.  (87)  as 


2n-l 

F2(«)QsvSV  ^zn) 


Substitution  of  Eq.  (84)  into  Eq.  (94)  gives 


2n-l 

Fi(w)F2(w)Qsvsv  ^i(^yn>  ^zn^  ^ 


Introducing  the  SV  wave  attenuation  constant  a  of  the  unidirec¬ 
tional  fiberglass  epoxy  composite  and  a  possible  electrical  sig¬ 
nal  amplification  factor  K,  Eq.  (95)  can  be  written  as 


V'  =  KF]^(w)F2(w)QgYgy  ^iC^yn*  ^zn^^ 


Eq.  (96)  gives  the  output  voltage  amplitude  from  the  non-contact 
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receiving  transducer  due  to  an  input  voltage  amplitude  V  at  the 
non- contact  transmitting  transducer  when  the  SV  wave  path  has 
undergone  n  reflections  from  the  bottom  face  of  the  unidirec¬ 
tional  fiberglass  epoxy  composite  plate  specimen,  as  shown  in 
Fig.  9. 
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DISCUSSION  AND  CONCLUSION 


In  the  acous to -ultrasonic  input -output  characterization  of  the 
unidirectional  fiberglass  epoxy  composite  plate,  the  angle  of 
reflection  of  the  reflected  SV  wave  is  equal  to  the  angle  of 
incidence  of  the  incident  SV  wave  for  each  reflection  at  either 
the  top  or  the  bottom  face  of  the  plate.  This  is  due  to  the  fact 
that  the  isotropic  plane  is  parallel  to  both  faces  of  the  plate. 
However,  if  the  parallelism  between  the  isotropic  plane  and  the 
plane  boundaries  where  reflection  occurs  does  not  exist,  the 
angle  of  reflection  is  not  equal  to  the  angle  of  incidence.  In 
such  a  case,  the  use  of  a  semi -infinite  transversely  isotropic 
medium,  neglecting  the  existence  of  the  bottom  face  of  the  plate 
except  for  the  cumulative  reflection  coefficients,  to  compute  the 
delay  time,  the  displacements  and  the  directivity  functions 
becomes  inappropriate . 

For  an  SV  wave  travelling  in  an  isotropic  medium,  its  phase 
velocity  is  a  constant  value;  that  is,  its  phase  velocity  is 
directionally  independent.  However,  as  shown  in  Fig.  10,  the 
phase  velocity  C;^  wave  travelling  in  a  unidirectional 

fiberglass  epoxy  composite  plate  is  a  function  of  the  angle  of 
incidence  6,  Consequently,  the  directional  dependence  of  the 
phase  velocity  of  an  SV  wave  in  the  unidirectional  fiberglass 
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epoxy  composite  plate  has  an  effect  on  the  delay  time  when  con¬ 
ducting  acousto-ultrasonic  input-output  characterization.  This 
phenomenon  can  be  observed  at  low  values  of  dimensionless  separa¬ 
tion  i/h  in  Fig.  11,  where  sharp  increases  in  delay  time  t^j  are 
caused  by  an  increase  in  the  number  of  reflection  n  from  the  bot¬ 
tom  face  of  the  plate  at  a  given  value  of  the  dimensionless  sepa¬ 
ration  i/h.  Low  values  of  dimensionless  separation  i/h  correspond 
to  an  angle  of  incidence  9  (given  in  Eq.  (49))  whose  value  is 
less  than  60°  which  is  the  abscissa  of  the  maximum  phase  velocity 
Cl,  as  shown  in  Fig.  10.  For  angles  of  incidence  9  less  than 
60°,  the  phase  velocity  increases  with  increasing  angle  of 
incidence  9,  as  shown  in  Fig.  10.  However,  an  increase  in  the 
number  of  reflections  n  from  the  bottom  face  of  the  plate  for  a 
given  plate  thickness  h  and  a  given  separation  distance  i  results 
in  a  decrease  in  the  angle  of  incidence  9  given  by  Eq.  (49), 
thereby  increasing  the  travelling  distance  according  to  Eq. 

(50)  and  decreasing  the  accompanying  phase  velocity.  The 
increase  in  the  travelling  distance  and  the  decrease  in  the 
phase  velocity  Cj^  account  for  the  sharp  increases  in  the  delay 
time  tj^  at  low  values  of  dimensionless  separation  t/h.  For  the 
case  of  an  SV  wave  travelling  in  an  isotropic  medium,  the 
increase  in  the  delay  time  tjj  caused  by  an  increase  in  the  number 
^®fl®ctions  n  is  attributed  solely  to  the  increase  in  the 
travelling  distance  Rj^. 
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Consider  the  case  of  a  given  number  of  reflections  n.  For  a 
given  plate  thickness  h,  an  increase  in  the  separation  distance  H 
results  in  an  increases  in  the  travelling  distance  R^,  as  shown 
in  Fig.  9,  and  an  increase  in  the  angle  of  incidence  9  given  by 
Eq.  (49).  An  increase  in  the  angle  of  incidence  9  may  result  in 
an  increase  or  a  decrease  in  the  SV-wave  phase  velocity 
depending  on  whether  the  angle  of  incidence  9  is  less  than  or 
greater  than  60\  as  shown  in  Fig.  10.  However,  as  the  separa¬ 
tion  distance  S,  approaches  infinity,  the  angle  of  incidence  9 
will  approach  90®,  and  the  phase  velocity  Ci  will  reach  a  limit, 
as  shown  in  Fig.  10.  Consequently,  (for  i/h  ®)  the  increase  in 
the  delay  time  t^  is  due  almost  exclusively  to  the  increase  in 
the  travelling  distance  R^,  which  is  similar  to  the  isotropic 
medium  case. 

This  theoretical  investigation  provides  a  step  forward  in  the 
quantitative  understanding  of  acousto -ultrasonic  nondestructive 
evaluation  (NDE)  parameters  such  as  the  stress  wave  factor  (SWF) 
in  transversely  isotropic  media.  It  also  provides  the  potential 
for  assisting  in  the  development  of  more  efficient  and  more 
revealing  NDE  schemes  utilizing  wave  propagation. 
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Fig.  2  Coordinate  system  (x,y,z)  in  the  analysis  of  single 
reflection  problem  at  stress-free  plane  boundary  of 
semi-infinte  transversely  isotropic  medium  for  angle 
of  incidence  less  than  critical  angle;  z  =  0  is  plane 
boundary  where  single  reflection  occurs,  z <  0  is  free 
space,  and  x =  0  is  plane  of  incidence. 
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Fig.  3  Two  sheets  of  slowness  surface  for  P  and  SV  waves  in  unidirectional  fiberglass 
epoxy  composite  for  positive  y— z  quadrant. 
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Fig.  4  Angle  of  reflection  of  reflected  P  wave  versus  angle  of  Incidence  of 
incident  SV  wave  in  unidirectional  fiberglass  epoxy  composite. 


ratio  of  reflected  SVwave 


Fig.  5  Amplitude  ratios  of  reflected  SV  and  P  waves  to  incident  SV  wav 
angle  of  incidence  in  unidirectional  fiberglass  epoxy  composite 
incidence  less  than  critical  angle. 


P  wave 

amplitude  decreases 
exponentially  with  z 


Fig.  6  Existence  of  reflected  surface  P  wave  whose  amplitude 

decreases  exponentially  with  z  for  angle  of  incidence  of 
incident  SV  wave  in  unidirectional  fiberglass  epoxy 
composite,  equal  to  or  greater  than  critical  angle  0^,; 
z =  0  is  plane  boundary  where  single  reflection  occurs, 
z <  0  is  free  space,  and  x =  0  is  plane  of  incidence. 
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:  z-component  energy  flux  of  ref  lected  SV  wove 

*  Z-componenf  energy  flux  of  reflected  Pwove 

*  Z-component  energy  ftux  of  incident  SV  wove 


Angle  of  incidence,  6i  (degree) 

ormal  to  plane  boundary,  assuming  radian  frequency  equal  to  unity, 
wave  in  unidirectional  fiberglass  epoxy  composite  for  angle  of 
than  critical  angle. 


=  500 
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bottom 


Schematic  illustrating  sinusoidal  point  load  exciting 
infinite  transversely  isotropic  medium,  where  xy  plane 
is  isotropic  plane  in  cartesian  coordinate  system 
defined  by  0(x,y,z). 
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Fig.  13  Polar  diagram  for  directivity  function  Dy^  of  shear  stress  Xyz  at  transmitting 

origin  associated  with  SV  wave  in  unidirectional  fiberglass  epoxy  composite  plate 
specimen  at  frequency  of  0.75  MHz  for  positive  y-z  quadrant  due  to  applied  point 
force  acting  along  y-direction. 


Fig.  14  Polar  diagram  for  directivity  function  of  shear  stress  x  at  transmitting 

origin  associated  with  SV  wave  in  unidirectional  fiberglass  epoxy  composite  plate 
specimen  at  frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point 
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Fig.  18  Polar  diagram  for  directivity  function  of  normal  stress  at  transmitt 
origin  associated  with  SV  in  unidirectional  fiberglass  epoxy  composite  plate 

at  frequency  of  2.25  MHz  for  positive  y-z  quadrant  due  to  applied  point  fore 
along  y-direction. 
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Polar  diagram  for  directivity  function  D  of  shear  stress  at  transmitting 

origin  associated  with  SV  wave  in  unidirectional  fiberglass  epoxy  composite  plate 
specimen  at. frequency  of  0.75  MHz  for  positive  y-z  quadrant  due  to  applied  point 
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Polar  diagram  for  directivity  function  of  shear  stress  T  ^  at  transmitting 

origin  associated  with  SV  wave  in  unidirectional  fiberglass  ^poxy  composite  plate 
specimen  at  frequency  of  2.25  MHz  for  positive  y~z  quadrant  due  to  applied  point 
force  acting  along  z-direction. 
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Directivity  function  D2Z  Q'ong  y  axis 

Polar  diagram  for  directivity  function  D^z  of  normal  stress  Tzz  transmitting  origin 
associated  with  SV  wave  in  unidirectional  fiberglass  epoxy  composite  plate  specimen  at 
frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point  force  acting  along 
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Directivity  function  D^z  along  y  axis 
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Fig.  24  Polar  diagram  for  directivity  function  of  normal  stress  Tzz  transmitting 

origin  associated  with  SV  wave  in  unidirectional  fiberglass  epoxy  composite  plate 
specimen  at  frequency  of  2.25  MHz  for  positive  y-z  quadrant  due  to  applied  point 
force  acting  along  z~direction. 


